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多 孔 介 质 PIAS 流 过 程 入 Q 效应 的 数值 研究 
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摘 要 : 本 文 利用 随机 泡沫 数目 守恒 模型 对 均 质 多 孔 介 质 内 泡沫 液 渗流 过 程 的 入 口 效应 进行 了 数值 研究 。 随 机 泡沫 数目 守恒 模 
型 仅 利用 两 个 参数 : 泡沫 生成 速率 K, 与 最 大 泡沫 个 数 rns 来 描述 多 孔 介 质 内 非 稳 态 驱 蔡 过 程 中 泡沫 的 生成 和 发 展 。 采 用 IMPES 
方法 对 二 维 介质 内 泡沫 两 相 渗流 过 程 的 控制 方程 组 进行 求解 ， 基 于 驱 替 过 程 中 液 相 压力 、 液 相 饱 和 度 及 气泡 密度 分 布 等 参数 的 
分 析 ， 对 泡沫 驱 蔡 过 程 的 入 口 效 应 进行 了 详细 研究 。 结 果 发 现 ， 增 大 泡沫 生成 速率 K,， 泡 沫 渗流 过 程 的 入 口 效 应 减 小 ; 而 增 
大 参数 xm， 虽 可 提高 驱 替 压 差 ， 但 对 多 孔 介质 入 口 段 各 参数 的 分 布 规律 影响 不 大 。 
关键 词 : 随机 泡沫 数目 守恒 模型 ， 多 孔 介 质 ;， 泡 沫 渗流 ;入口 效应 ;数值 研究 
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Abstract: Numerical investigation has been carried out concerning the inlet effect of foam flow process in a 
homogeneous porous media based on the stochastic bubble population balance model. The stochastic bubble 
population balance model could predict simply and robustly the bubble growth and development behavior in the 
transient foam displacement process in porous media by employing only two parameters: the foam generation rate of 
K, and the maximum bubble number of ns. The governing equations for two-dimensional two-phase foam seepage 
process are solved through IMPES method, and the distribution characteristics of water phase pressure, water 
saturation and bubble density are analyzed to investigate the inlet effect in a foam flow process in the porous media. 
The results show the parameter variations in the inlet region could be suppressed at higher K, values, indicating less 
inlet effect at higher foam generation rates. Although it shows higher pressure drop along the sample, increasing the 
nw Value has little effect on the distribution laws of various parameters in the entrance region of the porous media. 
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在 石油 开采 过 程 中 , 泡沫 技术 被 广泛 用 于 强化 
采油 (EOR) 的 工业 实践 。 传 统 的 气体 驱 油 过 程 由 于 
重力 分 层 及 气体 在 高 渗透 率 岩 层 的 窜 流 而 导致 较 
低 的 驱 油 效率 , 而 采用 泡沫 技术 则 可 显著 提高 气相 
黏度 ， 对 驱 蔡 气体 进行 良好 的 流 度 控制 ， 进 而 提高 
最 终 原 油 采 收 率 上 9。 

利用 数值 方法 可 简单 可 靠 地 对 多 孔 介 质 内 泡 
沫 的 驱 蔡 特性 进行 分 析 、 模 拟 和 预测 。 目 前 部 分 看 
究 者 采用 单 相 非 牛顿 流体 的 假设 来 模拟 泡沫 的 流 
ARE PEO, (AFTER ZS IAA PE, AS ABR ELIZ 
应 泡沫 在 多 孔 介 质 内 的 渗流 机 理 。 要 成 功 描述 泡沫 
流体 在 多 孔 介 质 中 的 驱 替 特性 , 则 必须 对 介质 内 泡 
沫 的 生成 及 潭 灭 行为 进行 准确 建 模 。 在 基于 两 相 流 
体 分 相 假设 的 流 变 学 模型 中 ， 泡 沫 数目 ne 是 非常 
重要 的 参数 。 现 今 的 泡沫 数目 模型 中 ， 最 著名 的 是 
分 相 流动 模型 及 气泡 数量 守恒 模型 。 分 相 流动 模型 
由 Rossen ”1 提出 ， 由 于 模型 中 没有 精确 考虑 气 
泡 数量 随时 间 的 变化 , 因此 多 用 于 模拟 泡沫 的 稳 态 
渗流 特性 。 而 Patzek 及 Kovscek 等 上 5 六 提出 的 气 
泡 数量 守恒 模型 ， 虽 然 基 于 气泡 密度 的 动态 平衡 ， 
但 由 于 模型 内 参数 过 多 ， 而 在 实践 中 难以 应 用 。 因 
此 目前 较为 成 功 的 油 藏 多 相 流 模拟 商业 软件 ， 如 
Eclipse”, CMG 公司 的 STARSP9 在 处 理 泡沫 流体 
时 并 不 是 通过 计算 泡沫 的 表 观 粘度 来 求 得 其 流 度 ， 
而 是 直接 对 泡沫 相 的 流 度 进行 经 验 性 修正 。 

最 近 ，Zitha 和 Du 请- 的 开发 了 一 种 针对 泡沫 
在 多 孔 介 质 中 运动 的 新 的 随机 泡沫 数目 守恒 模型 。 
泡沫 数目 守恒 模型 基于 以 下 几 点 基本 假设 : (1) 泡 沫 
作为 复杂 流体 , 它 的 流 变 学 特性 由 Herschel-Bulkley 
模型 描述 ; (2) 泡 沫 的 流 变 学 特性 本 质 上 取决 于 泡沫 
的 密度 (每 单位 体积 多 孔 介 质 内 泡沫 的 数目 ); (3) 最 
后 ,在 微观 水 平 上 ,泡沫 的 产生 看 作 是 一 个 随机 的 
过 程 ， 泡 沫 产生 动力 学 遵循 简单 的 指数 增长 定律 。 
相 比 传统 的 数量 平衡 模型 ， 随 机 泡沫 数量 守恒 模型 
仅 包 括 两 个 基本 变量 ， 因 此 更 易于 由 实验 来 确定 。 

本 文采 用 随机 泡沫 数目 守恒 模型 对 泡沫 液 在 均 
质 多 孔 介质 内 渗流 过 程 的 入 口 效 应 进行 了 数值 模 
拟 。 所 谓 的 入 口 效应 ， 是 指 在 入 口 处 一 段 区 域内 ， 
某 些 变量 的 分 布 与 其 他 区 域内 的 这 些 变 量 的 分 布 
有 较 大 的 差异 。 在 泡沫 驱 替 过 程 中 ， 这 些 变量 包括 
液 相 饱和 度 ， 液 相 压力 以 及 泡沫 密度 。 本 文采 用 
IMPES 方法 对 二 维 两 相 泡沫 渗流 过 程 的 控制 方程 
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组 进行 求解 ， 基 于 驱 蔡 过 程 中 液 相 压力 、 液 相 饱 和 
度 及 气泡 密度 分 布 等 参数 的 分 析 , 对 泡沫 驱 替 过 程 
的 入 口 效应 进行 了 详细 研究 。 

1 计算 

1.1 计算 区 域 

图 1 为 本 文 计算 区 域 , 是 长 度 为 L， 宽 度 为 d， 
IREKE, 渗透 率 为 K 的 二 维 均 质 多 孔 介 质 。 流 
体 从 左边 中 心 部 分 流入 ， 从 右边 边界 流出 。 开 始 时 
先 将 介质 中 注 满 水 , 即 介质 中 液 相 饱 和 度 为 100%， 
然后 气体 和 表面 活化 剂 溶液 以 速度 u 从 入 口 边界 
IA, Husu, +u, U, Mu, 分 别 气体 和 液体 的 
分 速度 ， 泡 沫 在 介质 内 生成 、 发 展 并 对 其 中 原 有 液 
相 进 行 驱 替 ， 最 后 由 出 口 边界 流出 。 出 口 压力 保持 
不 变 为 0.1MPa。 
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图 1. 计算 区 域 

Fig.1. Calculation area 
1.2 基于 随机 泡沫 数目 守恒 模型 的 控制 方程 
泡沫 在 多 孔 介质 内 驱 奉 过 程 的 控制 方程 为 : 


ol aS, 
pS). 9 (pu,)=0 
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u, = 一 VB 
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S,+S,=; P=P.-P, (2) 
Po So U,> A, AP AAA i REE, 
饱和 度 , 速度 , 流 度 和 压力 , 其 中 ie {w, /1( 下 标 f 
指 的 是 泡沫 )。 控 制 方程 (1) 中 第 三 式 即 为 随机 泡沫 
数目 守恒 模型 ， 其 中 泡沫 生成 速率 Ks 和 最 大 泡沫 
HA n, 是 与 泡沫 产生 与 发 展 相关 的 两 个 参数 ， 可 
由 实验 测 得 。 计 算 过 程 中 假设 液 相 不 可 压缩 ， 而 泡 
沫 相符 合理 想 气体 状态 方程 。 辅 助 方程 中 p, AE 
细 压 力 ， 其 仅 为 液 相 饱 和 度 Sw 的 函数 ， 即 : 


1 1 
= p 0.5 = Swe A Sy, — Sre a (3) 
p= ar {SSS |S 


其 中 on , 指 的 是 气 液 两 相 之 间 的 表面 张力 ，” 


we 


为 拟 合 参 数 ， 为 第 数 。 


第 ; 相 流体 的 流 度 九 定义 为 ; 
ja (4) 
MM, 
其 中 为 第 i MATER, 为 流 
体 的 粘度 ， 泡 沫 流体 的 粘度 由 下 式 计算 
My = Hy + a— (5) 
: i: 
Rta, HIEM 
体积 多 孔 介质 内 泡沫 的 个 数 )。 
流体 的 相对 渗透 率 , 的 表达 式 为 : 


n 为 泡沫 密度 (每 单位 
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1.3 数值 求解 方法 

本 文采 用 IMPES 方法 3] 对 控制 方程 组 (1 
行 求解 。 此 问题 为 二 维 两 相 渗流 过 程 ， 数 值 求 解 的 
过 程 分 为 三 个 步骤 (上 标 n 代表 时 间 ): 

(1) 第 n 时 刻 的 压力 场 P? 由 前 一 时 刻 的 压力 场 
P, RR S 和 泡沫 密度 n” 隐 式 计算 得 出 ， 第 
n 时 刻 的 压力 场 P 可 以 由 P, AP" 之 和 求 得 其 


= 


H P" 由 方程 (3) 得 到 。 

(2) 第 ntl 时 刻 的 饱和 度 分 布 S%” 可 以 由 S% 和 
PP BARA. 

(3) 第 n+l 时 刻 的 泡沫 密度 分 布 a"”” 也 可 以 用 
Caen” | Pe S SH ERRE. 
1.4 数值 模拟 中 的 参数 
数值 模拟 中 采用 的 参数 如 表 1 和 表 2 所 示 。 

表 1 数值 模拟 参数 


Table 1 Numerical simulation parameters 


Ski “(sis 参数 ë 数值 
KÆ Lm) 0.05 Ug in( 10/8) 0.0002 
宽度 d(m) 0.025 Wwin(m/s) 0.00002 
时 间 间 隔 和 Az(s) 0.0001 Vain(m/s) 0 
义 方向 网 格 数 N 52 Vwin(m/S) 0 
立方 向 网 格 数 M 12 PoPa) 101325 
时 间 z(s) 500 


表 2 流体 与 岩石 特 ' 性 参数 


Table 2 Fluid and rock characteristics parameters 


u,(Pa*s) 1004.0 10° 
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参数 数值 参数 数值 
Sue 0.1 (Pa ` s) 18.1x10° 
Q 0.23 a 8.6X107 
Km 2.0X102 A 5.0 
ky 1.0 c 0.333333 
ki, 0.5 Y 0.5 


0 w( N/m) 30 10° 
YE: 表 2 中 参数 参考 Bentheimer At»), 


2 数值 模拟 结果 及 分 析 


2.1 泡沫 生成 速率 K, 的 影响 (7-=200) 

2.1.1 参数 二 维 分 布 云图 

图 2、 图 3 分 别 为 ns=200 保持 不 变 、 驱 蔡 过 
程 完成 后 同一 时 刻 (T=500s)、 不 同 K, 值 时 液 相 饱 和 
度 、 液 相 压力 和 泡沫 密度 分 布 的 二 维 云图 。 从 图 中 
可 以 看 出 ,在 多 孔 介质 入 口 部 分 ， 液 相 饱 和 度 在 ， 
心 区 域 的 数值 要 大 于 两 侧 , 介质 中 心 的 液 相 压 力也 
明显 大 于 两 侧 ， 泡 沫 密度 在 入 口 处 最 小 , 但 沿 程 逐 
渐 增 加 到 最 大 泡沫 数目 。 入 口 段 参数 分 布 的 明显 不 
同 表 明 泡 沫 液 在 多 孔 介 质 驱 蔡 过程 中 具有 明显 的 
入 口 效应 。 比 较 图 2 与 图 3 还 可 看 出 ， 图 2 在 介质 
Is 数 的 变化 范围 明 
大 于 图 3， 说 明 KERK, 其 入 口 效应 越 小 。 这 
er ais K, 越 大 ， 说 明 泡 
沫 生成 越 快 并 能 够 在 较 小 的 范围 内 达到 充分 发 展 ， 
此 其 驱 蔡 过 程 的 入 口 效应 较 小 。 


a engin 
2 K=0.1 MERKER ESA h RAEE TF AR 
和 度 、 液 相 压 力 及 泡沫 密度 分 布 图 ) 
Figure 2 Distribution diagram of various displacement 
parameters at K,=0.1 (From top to bottom are water saturation, 
water pressure and foam density distributions) 
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3 K=1.0 时 泡沫 驱 蔡 特性 参数 分 布 图 (从 上 至 下 为 液 相 饱 
和 度 、 液 相 压 力 及 泡沫 密度 分 布 图 ) 
Figure 3 Distribution diagram of various displacement 


parameters at K,=1.0 (From top to bottom are water saturation, 
water pressure and foam density distributions) 


2.1.2 驱 蔡 过 程 参 数 沿 程 分 布 特性 

图 4、 图 5 分 别 为 动态 驱 蔡 过 程 中 不 同 K。 值 时 
液 相 饱 和 度 、 液 相 压力 和 泡沫 密度 沿 介 质 长 度 方向 
分 布 曲线 。 沿 程 参数 由 同一 横 截 面 上 的 参数 平均 获 
得 。 图 中 可 清楚 看 到 两 种 K, 值 下 泡沫 液 在 200s 时 
均 已 在 介质 出 口 达到 突破 。 由 图 中 还 可 看 出 ， 两 种 
K, IE F, 泡沫 完成 突破 后 介质 入 口 区 域 液 相 饱 和 度 
数值 较 大 ， 液 相 压力 降低 幅度 较 小 ， 泡 沫 数目 变化 
明显 ， 均 显示 出 明显 的 入 口 效应 。 对 比 图 4 与 图 5 
可 以 看 出 ， 图 4 的 入 口 效应 更 加 明显 。 说 明 K, 
大 ， 泡 沫 在 多 孔 介 质 入 口 区 域 的 参数 变化 范围 越 
小 ， 即 入 口 效应 越 小 ， 其 结论 与 二 维 云图 分 析 结果 
一 致 。 
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Figure 4 Distributions of various foam displacement 
parameters along the flow direction at K,=0.1 (From top to 
bottom are water saturation, water phase pressure and foam 


density) 
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Figure 5 Distributions of various foam displacement 

parameters along the flow direction at K,=1.0 (From top to 

bottom are water saturation, water phase pressure and foam 


density) 

2.2 最 大 泡沫 个 数 me 的 影响 (K=0. 1) 

2.2.1 参数 二 维 分 布 云图 

图 6 驱 替 过 程 完成 后 (T=500s)、mx=500 时 液 相 
饱和 度 、 液 相 压 力 和 泡沫 密度 分 布 的 二 维 云 图 。 从 
图 中 可 以 看 出 在 多 孔 介 质 入 口 区 域 , 试 样 中 心 部 分 
的 液 相 饱 和 度 要 大 于 两 侧 ， 液 相 压力 也 大 于 两 侧 ， 
泡沫 密度 在 入 口 处 最 小 。 而 在 试 样 后 半 部 分 三 个 参 
数 的 变化 均 明 显 减 小 , 因此 驱 替 过 程 也 显示 出 明显 
的 入 口 效 应 。 与 图 2 n200 的 结果 对 比 ， 虽 然 图 6 
在 xx 较 大 时 具有 更 小 的 液 相 饱 和 度 、 更 大 的 驱 替 压 
差 及 更 多 的 泡沫 数目 , 但 两 图 中 泡沫 驱 蔡 过 程 参 数 


PS 


空间 分 布 规律 基本 相同 , 说 明 nw 的 大 小 对 泡沫 液 在 
多 孔 介 质 内 驱 蔡 过 程 的 入 口 效 应 影响 不 大 。 
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6 nw=500 时 泡沫 驱 蔡 特性 参数 分 布 图 (从 上 至 下 分 别 为 液 


Ds 


相 饱 和 度 、 液 相 压 力 及 泡沫 密度 ) 
Figure 6 Distribution diagram of various displacement 
parameters at n„=500 (From top to bottom are water saturation, 
water pressure and foam density distribution) 


2.2.2 驱 替 过 程 参 数 沿 程 分 布 特性 

图 7 为 we=500 时 液 相 饱和 度 、 液 相 压 力 和 泡 
沫 密度 沿 介质 长 度 方向 的 分 布 曲线 。 由 图 中 可 以 看 
出 ，nw=500 时 多 孔 介 质 前 部 的 液 相 饱和 度 具 有 下 
降 趋势 、 多 孔 介 质 入 口 区 域 的 压 降幅 度 较 介质 内 后 
部 区 域 小 、 介 质 入 口 段 的 泡沫 密度 变化 更 为 显著 。 
因此 mx=500 时 泡沫 驱 蔡 过 程 也 具有 明显 的 入 口 效 
应 。 进 一 步 与 图 4 n=200 的 结果 对 比 可 以 看 出 ， 
两 种 情况 下 虽然 各 参数 的 绝对 数值 并 不 相同 , 但 各 
参数 的 沿 程 分 布 规律 基本 一 致 , 说明 me 取 值 对 泡沫 
液 在 多 孔 介 质 内 驱 蔡 过 程 的 入 口 效 应 影响 并 不 明 
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7 ns=500 时 不 同时 刻 泡 沫 驱 替 特性 参数 沿 流动 方向 的 分 
布 (从 上 至 下 为 液 相 饱 和 度 、 液 相 压力 及 泡沫 密度 ) 
Figure 7 Distributions of various foam displacement 
parameters along the flow direction at n,=500 (From top to 
bottom are water saturation, water phase pressure and foam 


density) 


3 结论 


本 文采 用 随机 泡沫 数目 守恒 模型 SNA 
二 维 均 质 多 和 孔 介 质 内 渗流 过 程 的 入 口 效 应 进行 了 
数值 模拟 分 析 ， 得 出 如 下 结论 : 

d) 泡沫 由 中 心 注 入 多 孔 介 质 后 ， 在 介质 内 生 
成 并 逐步 发 展 ， 其 驱 替 过 程 显示 出 明显 的 
入 口 效应 ; 

(2) ”泡沫 生成 速率 参数 有 增 大 ， 泡 沫 液 在 多 和 孔 
介质 入 口 区 域 的 参数 变化 范围 变 小 ， 即 入 
口 效应 减 小 ; 

G) ”最 大 泡沫 密度 参数 zw 不 同 虽 可 影响 各 驱 替 
参数 的 绝对 数值 ， 但 各 参数 的 沿 程 分 布 规 

律 基本 一 致 ， 说 明 me 取 值 对 泡沫 驱 替 过 程 

入 口 效 应 的 影响 并 不 明显 。 


也 
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